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ABSTRACT 

A theoretical study of Ali_ x Ga x N, based on full-potential linearized augmented plane wave 
method, is used to investigate the variations in the bandgap, optical properties and non-linear 
behavior of the compound with the variation of Ga concentration. It is found that the bandgap 
decreases with the increase of Ga in Ali_ x Ga x N. A maximum value of 5.5 eV is determined for 
the bandgap of pure A1N which reaches to minimum value of 3.0 eV when Al is completely 
replaced by Ga. The static index of refraction and dielectric constant decreases with the increase 
in bandgap of the material, assigning a high index of refraction to pure GaN when compared to 
pure A1N. The refractive index drops below 1 for photon energies larger than 14 eV results group 
velocity of the incident radiation higher than the vacuum velocity of light. This astonishing result 
shows that at higher energies the optical properties of the material shifts from linear to non- 
linear. Furthermore, frequency dependent reflectivity and absorption coefficients show that peak 
value of the absorption coefficient and reflectivity shifts towards lower energy in the UV 
spectrum with the increase in Ga concentration. This comprehensive theoretical study of the 
optoelectronic properties of the alloys is presented for the first time which predicts that the 
material can be effectively used in the optical devices working in the visible and UV spectrum. 
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I. INTRODUCTION 

Group-Ill nitride family and their alloys are suitable materials for optoelectronic 
applications in short wave length, high temperature, high frequency and high power electronic 
devices [1,2]. Rare-earth doped A1N and GaN are effectively used in photonic devices, 
electroluminescent devices (ELDs), diode lasers, medicine and health physics [3-9]. Accurate 
knowledge of the bandgap and optical properties of these compounds is very important for the 
design and analysis of various optoelectronic and photonic devices. It is important to understand 
fundamental optical properties over a wide range of wavelengths. In reality, optical properties 
reflect the density of states (DOS) of a compound and their analysis is one of the most effective 
tools to understand the electronic structure of the material [10, 11]. 

Ali_ x Ga x N and Ali_ x In x N are well suited materials as photo-detectors in the UV spectrum 
[12]. Many experimental and theoretical studies have been performed on the wurtzite structure of 
these ternary alloys [13-15], but little work has been reported on the zinc -blend phase of these 
materials [16], even after the successful growth reported by Okumura et al. [16] and Martinez et 
al. [17]. The investigation of zinc -blend structure of these materials is more important than 
wurtzite structure due to its large optical gain and lower threshold current density [18]. 

Lee and Wang [19] studied the electronic structure of zinc -blend Ali_ x Ga x N by local 
density approximation (LDA) and screened-exchange local density approximation (sX-LDA). 
While Kanoun et al. [20] reported the structural parameters and bandgap bowing in zinc -blend 
phase of Ali_ x Ga x N and Ali_ x In x N using density functional theory with local density 
approximation (LDA). It has been confirmed that generalized gradient approximation (GGA) is 
far more superior to LDA in the band structure calculation [21, 22], that's why GGA is used in 
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the present work. Though bandgap and electronic structure of Ali. x Ga x N has already been 
studied in Ref. [19] and Ref. [20] but no theoretical work is reported on the optoelectronic 
properties of the compound. Optical properties of Ali_ x Ga x N and their dependence on the 
variation of Ga concentration are calculated. 

II. THEORY AND CALCULATION 

The dielectric function s(a>) of the electron gas, with its strong dependence on frequency 
has significant effect on the physical properties of solids. It describes the collective excitations of 
the Fermi Sea, such as the volume and surface plasmons. The dielectric function depends on the 
electronic band structure of a crystal and its investigation by optical spectroscopy is a powerful 
tool in the determination of the overall band behavior of a crystal. It can be divided into two 
parts; real and imaginary parts: 

s(co) = e x (a>) + is 2 (a>) ( 1 ) 

The imaginary part of the complex dielectric function s 2 (&>) in cubic symmetry compounds can 
be calculated by the following relation [23, 24]: 

- 2 (-)=A^Ji p -w 2 ^\i (2) 

27KO ~ i z VcO nn \k) 

s 2 (co) is strongly dependent on the joint density of states, co nn , , and momentum matrix element, 
P nn , . While Kramers-Kronig relation is used to calculate the real part of the dielectric function 
s^co) from £ 2 (co) [25]: 




The values of the real and imaginary parts of the frequency dependent dielectric function 
provides basis for the calculation of the refractive indices n (a>) . The refractive indices of ternary 
alloy are important in the designing of the optical parameter [26]. The complex refractive index 
of a compound is: 

n(a>) = n(co) + ik(co) = s V2 = (s l +is 2 ) U2 (4) 
Where n(co) represent the real part of the refractive index and k(a>) represent the imaginary part 
or the extinction coefficient [27]. The real part of the refractive index n(co) and the extinction 
coefficient k(co) can be calculated from and s 2 (co) using the following equations: 
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Using the above optical parameters n(co) and k(co) ; the frequency dependent reflectivity, R{co) 
can be calculated by the following relation: 



R(co) = 



n-l 



n+l 
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Similarly the absorption coefficient is given by Beer's law [27]: 

2kco 47ik 



a = 



(8) 



The absorption coefficient can also be calculated from the dielectric function [27]: 

_ i/ 

a{co) = 2cok(a>) = 42co {s la>) 2 + s 2 (co) 2 }^ 2 - s x (co) 



(9) 



In this article, all density functional calculations are performed with the full-potential 
linearized augmented plane wave method, using Wu-Cohen potential in the generalized gradient 



approximation embedded in the wien2k software [28]. The details of calculation are reported in 
our previous work [5, 6, 24, 24, 29]. In the present work we used Rmt * K max = 8.00 value for 
the plane wave cut-off in the interstitial region. For the calculation of the optical properties of 
Ali_ x Ga x N; a mesh of 3500 k-points is taken for the Brillouin zone integrations in the 
corresponding irreducible wedge. 
III. RESULTS AND DISCUSSION 

Electronic bandgaps, optical parameters such as dielectric constants and refractive indices 
for Ali_ x Ga x N are calculated. Fig.l shows that the calculated bandgap varies linearly as a 
function of Ga concentration. Bandgap is of direct nature for the full compositional range of Ga 
in Ali_ x Ga x N (0 < x < 1). For x = with no Ga concentration, the bandgap of pure A1N is 
determined to be 5.50 eV. By increasing the fractional concentration x of Ga, the bandgap of the 
materials decreases. This decrease saturates at the maximum value of x and gives a bandgap of 
3.00 eV for pure GaN. Moreover, our calculated results of the bandgaps for Ali_ x Ga x N (0 < x < 
1) is shown in Fig.l, are closer to the experimental results of A1N and GaN [30, 31]. It is further 
noted that our theoretical results are much improved than reference [20]. The reason for our 
better results is using Wu-Cohen potential in the GGA scheme [32] which is more effective in 
the bandgap calculation than the reported results of LDA [21, 22]. 

The variation in the bandgap with the constituents' concentration occurs due to the 
difference in the density of states [6]. It has already been reported that the density of states 
affects the band structure of a compound and hence its bandgap [33-35]. Furthermore, difference 
in the densities of states of A1N and GaN are responsible for making A1N as a wider bandgap 
material than GaN. When some of the Al atoms are replaced by Ga atoms in the A1N crystal, 
then the mixing of Al and Ga in the compound is responsible for the change in the density of 
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states. At a higher Al concentration the distribution of the density of states of Al is dominant, 
making the bandgap of the resulting material closer to the bandgap of pure A1N. In a similar way 
at higher Ga concentration, the bandgap of Ali_ x Ga x N is closer to pure GaN due to the dominant 
effect of the distribution of the density of states of GaN. 

The variation in the bandgap of Ali_ x Ga x N provides a promising result of the bandgap 
engineering for optoelectronic devices. It covers a full spectrum of wavelengths from visible to 
UV, because the bandgap of pure A1N is transparent to visible light spectrum and part of the 
ultraviolet light while pure GaN is transparent to visible light. Depending upon the need and 
requirement of a particular application, any desired bandgap in this material can be achieved 
within a limit of 3.0 eV to 5.50 eV. Thus by varying the Ga concentration in Ali_ x Ga x N, 
optoelectronic devices in various regions of the spectrum can be formed. It can also be used to 
construct quantum wells of a desired width and depth [6, 36, 37]. 

The calculated real and imaginary parts of the dielectric function of Ali_ x Ga x N at x = 0, 
0.25, 0.50, 0.75 and 1.0 for the energy range to 25 eV are shown in Fig.2. It is clear from the 
figure that for x = 0, 0.25, 0.50, 0.75 and 1.0 the critical points in the imaginary parts of the 
dielectric function occurs at about 5.43 eV, 4.70 eV, 4.19 eV, 3.53 eV and 2.86 eV. These points 
are closely related to the direct band gaps (E^~ T ); 5.50 eV, 4.78 eV, 4.30 eV, 3.60 eV and 3.0 
eV of Ali_ x Ga x N at x = 0, 0.25, 0.50, 0.75 and 1.0. It is clear from the figure that A1N have strong 
absorption region 7.01 eV to 15.35 eV. This region consists of different peaks. The width of the 
absorption region decreases by changing the concentration of Ga from to 100 % and the 
sharpness and height of the peaks increases. The variations in the linear optical absorption region 
from one compound to another can be related to the bandgap of the materials. The materials with 
a bandgap lesser than 3.1 eV work well in the visible light devices applications while those with 
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bandgap larger than 3.1 eV can also be used in ultraviolet (UV) applications [6-9, 5]. These 
prominent variations in the optical absorption region with bandgap of Ali_ x Ga x N (3.0 eV to 5.50 
eV) confirms its suitability for device applications in the major parts of the spectrum; visible to 
UV. 

Fig. 2 also show the calculated real parts of the complex dielectric function s^co) for Ali_ 
x Ga x N. It is clear from the figure that the static dielectric constant s x (o) , the low energy limit of 
s { (eo), is strongly dependent on the bandgap of a compound. The calculated values of £ - 1 (o)Vs 
bandgap is plotted in Fig. 3. The plot shows inverse relation of £ - 1 (o)with the bandgap. The 
inverse relation of £ - 1 (o)with the bandgap can be explained by Penn Model [38]. 

£(0)*l + (/^/E g ) 2 (1Q) 

The above relation can be used to calculate E g using the value of £ - 1 (o)and plasma energy %co p . It 

is evident from Fig.2 that the first hump for A1N appears at 7.16 eV. The hump is shifted toward 
higher energies by increasing the concentration of Ga. 

For all the investigated compounds ofAli_ x Ga x N at x = 0, 0.25, 0.50, 0.75 and 1.0; s^ 00 ) 
becomes zero at certain energy and then decreases to minimum values, with negative numbers, 
about 14.98 eV, 14.35 eV, 14.21 eV, 13.58 eV and 13.45 eV. The negative values of ^(<y)show 
that in this energy region the incident electromagnetic waves are totally reflected from the 
medium, hence the material exhibit metallic nature. The negative values of s^co) correspond to 
the local maxima of the reflectivity (can be confirmed by Fig. 6). On further increase in the 
energy it increases and becomes zero at high energy limits. 

The calculated refractive index and extinction coefficient of Ali_ x Ga x N at x = 0, 0.25, 0.50, 0.75 
and 1.0 are presented in Fig. 4. A broad spectrum for n(a>) over a wide energy and frequency 
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range is noted for these compounds. The spectrum of n{co) closely follow s x (co) [27]. It is clear 
from the figure that the refractive index of the material increases with the increase in the Ga 
concentration in A1N. Fig. 5 shows the plot between n(o)and bandgap. The plot confirms the 
inverse relation of refractive index with the bandgap of Ali_ x Ga x N. Two different features can be 
observed in the refractive index in Fig. 4. First; a maxima can be observed in the form of a bump 
in the spectrum at a particular energy and, second; the maxima shifts to lower energy region with 
the increase in Ga concentration. Pure A1N with x = has the broadest spectrum of n(co)with a 

maximum value around 2.0 at 8.0 eV. At intermediate energies a few bumps appear and then the 
curves vanish at higher energies. The reason for these vanishing curves at higher energy or 
frequency is due to the fact that, beyond certain energy, the material can no longer act as 
transparent material and it absorbs high energy photons. Perhaps the high bandgap of A1N makes 
it transparent up to high energy and frequency range [9]. However by increasing the 
concentration of Ga the appeared bumps in A1N changes to sharp peaks. The peaks are clear in 
pure GaN. At this point the refractive index reaches to its maximum value of more than 2.0 for 
the peak. Beyond this peak, the value of n(a>) dissipates rapidly. GaN absorbs most of the high 
and intermediate energy photons and is transparent in the low energy region due to the smaller 
bandgap of GaN than A1N [8,9]. 

It is also clear from Fig.4 that the refractive index falls below unity at certain frequencies. 
Refractive index lesser than unity (v g =c/n) shows that the group velocity of the incident radiation 
is greater than c. It means that the group velocity shifts to negative domain and the nature of the 
medium changes from linear to non-linear. In other words the material becomes superluminal for 
high energy photons [39, 40]. 
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A relationship between k{co) and energy is also shown in Fig. 4. The response of k(co) to 
the varying constituent's concentrations of the Ali_ x Ga x N is closely matching to the response of 
s 2 {co). This result shows a good agreement with the available literature [27]. The peak value 
shifts to lower energies as x increases from to 1 . 

Frequency dependent reflectivity, R(co), versus energy for Ali_ x Ga x N, is shown in Fig. 6. 
Peaks in the figure show that for each concentration there is a maximum value of reflectivity. 
The maximum lie in the energy range 13 - 20 eV and arises from the inter band transition. While, 
the minimum of reflectivity occurs in the energy range 8 - 10 eV is due to the collective plasma 
resonance. The depth of the plasma resonance can be determined by the imaginary part of the 
dielectric function [41]. A more interesting aspect of the figure is the shift of the peak value 
towards lower energies with the increase in the Ga concentration. The material possesses 
reflectivity in wide energy and frequency range. For pure A1N the maximum value of reflectivity 
is 0.5, occurs at 17 eV and reduces to 0.46 at 13 eV for pure GaN. Hence, the reflectivity of Ali_ 
x Ga x N varies with the Ga concentration. This characteristic of the compound makes it suitable 
for Bragg's reflector in various wave lengths [42, 43]. Since the peak values of reflectivity in Ali_ 
x Ga x N changes with Ga concentrations, hence Bragg's reflector in the desired wavelength can be 
made and controlled by the Ga concentrations. 

Frequency dependent absorption coefficient for Ali_ x Ga x N (0 < x < 1) is shown in Fig.7. 
The variation in frequency is expressed in terms of energy; ranging from up to 30 eV. From the 
figure it is clear that the peak correspond to maximum absorption coefficient is shifted towards 
lower energies as the Ga concentration increases. Furthermore the maximum value drops from 
3 00* 10 4 cm" 1 (for pure A1N) to 250*10 4 cm _1 (for pure GaN). From the figure, two interesting 
features of the absorption coefficient can be noted; one is the maximum value of the absorption 
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coefficient shifts from 16 eV in pure A1N to 13 eV in pure GaN, while the second interesting 
feature is the appearance of the second peak with the increase in Ga concentration. The second 
peak is not observed in pure A1N, while two peaks are clear in GaN. These two peaks in the 
absorption coefficient explain the absorption of maximum light at two different wave lengths. 
Due to this property of the material, it can be used for wavelength filtering purposes in those 
regions. 

IV. CONCLUSION 

Density functional calculations are used to study optoelectronic properties of Ali. x Ga x N 
(0 < x < 1). Optical parameters like index of refraction, dielectric constant, reflectivity and 
absorption coefficient are strongly dependent on the concentration of Ga in the Ali. x Ga x N 
crystal. Zero frequency limit of dielectric function, refractive index and reflectivity decreases 
with the increase of Ga in Ali_ x Ga x N. The material becomes superluminal for photon energies 
larger than 14 eV, where the refractive index drops below 1. The peak value of the absorption 
coefficient and reflectivity shifts towards lower energy with the increase in Ga concentration, 
while two prominent peaks appear in the imaginary part of the dielectric function with Ga in A1N 
and the peak values enhances with the increase in Ga. The prominent variations in the optical 
parameters in the energy range 3 - 15 eV with Ga concentration in Ali_ x Ga x N makes it suitable 
for optical devices in the major parts of visible and UV spectrum. 
Acknowledgements 

Prof. Dr. Keith Prisbrey, MSE, University of Idaho and Prof. Dr. Nazma Ikram, Ex. 
Director, Center for Solid State Physics, Punjab University are highly acknowledged for their 
valuable suggestions. 



11 

REFERENCES 

[I] M. Esmaeili, M. Gholam, H. Haratizadeh, B. Monemar, P. O. Holtz, S. Kamiyama, H. 
Amano, and I. Akasaki, Opto-Electronics Review 17, 293 (2009). 

[2] H. X. Jiang, and J. Y. Lin, Opto-Electronics Review 10, 271 (2002). 
[3] F. A. Ponce, and D. P. Bour, Nature 386, 351 (1997). 

[4] R. Hui, S. Taherion, Y. Wan, J. Li, S. X. Jin, J. Y. Lin, and H. X. Jiang, Appl. Phys. Lett. 
82, 1326 (2003). 

[5] B. Amin, I. Ahmad, and M. Maqbool, Journal of Light Wave Tech. 28, 223 (2010). 
[6] M. Maqbool, B. Amin, and I. Ahmad, JOSA B, 26, 2180 (2009). 
[7] M. Maqbool, and I. Ahmad, Current Applied Physics 9, 234 (2009). 
[8] M. Maqbool, M. E. Kordesch, and I. Ahmad, Current Applied Physics 9, 417 (2009). 
[9] M. Maqbool, I. Ahmad, H. H. Richardson, and M. E. Kordesch, Appl. Phys. Lett. 91, 
193511 (2007). 

[10] J. Wu, E. E. Haller, H. Lu, W. J. Schaff, Y. Saito, and Y. Nanishi, Appl. Phys. Lett. 80, 
3967 (2002). 

[II] T. Matsuoka, H. Okamoto, M. Nakao, H. Harima, and E. Kurimoto, Appl. Phys. Lett. 81, 
1246 (2002). 

[12] T. Li, J. C. Carrano, C. J. Eiting, P. A. Grudowski, D. J. H. Lambert, H. K. Kwon, R. D. 

Dupuis, J. C. Campbell, and R. T. Tober, Fiber and Integrated Optics 20, 125 (2001). 
[13] S. Nakamura, "The Blue Laser Diode GaN Based Light Emitters and Lasers" (Springer, 

Berlin, 1997). 

[14] Z. Dridi, B. Bouhafs, and P. Ruterana, New Journal of Physics 4, 94.1 (2002). 



12 

[15] S. R. Lee, A. F. Wright, M. H. Crawford, G. A. Petersen, J. Han, and R. M. Biefeld, 

Appl. Phys. Lett. 74, 3344 (1999). 
[16] H. Okumura, H. Hamaguchi, T. Koizumi, K. Balakrishnan, Y. Ishida, M. Arita, S. 

Chichibu, H. Nakanishi, T. Nagatomo, and S. Yoshida, J. Cryst. Growth 189, 390 (1998). 
[17] E. Martinez-Guerrero, E. Bellet-Amalric, L. Martinet, G. Feuillet, B. Daudin, H. 

Mariette, P. Holliger, C. Dubois, C. Bru-Chevallier, P. Nze, T. Aboughe Chassagne, 

G. Ferro, and Y. Monteil, J. Appl. Phys. 91, 4983 (2002). 
[18] S. H. Park, and S. L. Chuang, J. Appl. Phys. 87, 353 (2000). 
[19] B. Lee and L. W. Wang, Phys. Rev. B 73, 153309 (2006). 

[20] M. B. Kanoun, S. Goumri-Said, A. E. Merad, and H. Mariette, J. Appl. Phys. 98, 063710 
(2005). 

[21] P. Dufek, P. Blaha, V. Sliwko, and K. Schwarz, Phys. Rev. B 49, 10170 (1994). 
[22] P. Dufek, P. Blaha, and K. Schwarz, Phys. Rev. B 50, 7279 (1994). 
[23] M. A. Khan, A. Kashyap, A. K. Solanki,T. Nautiyal, and S. Auluck, Phys. Rev. B 48, 
16974 (1993). 

[24] B. Amin, I. Ahmad, M. Maqbool, N. Ikram, Y. Saeed, A. Ahmad, and S. Arif, Journal of 

Alloys and Compounds, 70, 874 (2010). 
[25] F. Wooten, "Optical Properties of Solids," Academic Press, New York, 1972. 
[26] M. J. Bergmann, U. Ozgur, H. C. Casey, H. O. Everitt, and J. F. Muth, Appl. Phys. Lett. 

75, 67 (1999). 

[27] M. Fox, "Optical Properties of Solids," Oxford University Press, 2001 . 



13 

[28] P. Blaha, K. Schwarz, G. K. H. Madsen, D. K.vanicka, and J. Luitz, "WIEN2K-An 

Augmented Plane Wave and Local Orbital Program for calculating Crystal Properties. 
Techn. Universitate Wien, Austria," ISBN:3-9501031-l-l-2. (2001). 

[29] B. Amin, and Iftikhar Ahmad, J. Appl. Phys. 106, 0937108 (2009). 

[30] W. J. Fan, M. F. Li, and T. C. Chong, J. Appl. Phys. 79, 188 (1996). 

[31] S. Bloom, G. Harbeke, E. Meier, and I. B. Ortenburger, Phys. Status. Solidi, B 66, 161 
(1974). 

[32] Z. Wu, and R. E. Cohen, Phys. Rev. B. 37, 2351 16 (2006). 

[33] L. C. Duda, C. B. Stagarescu, J. Downes, K. E. Smith, D. Korakakis, T. D. Moustakus, J. 

Guo, and J. Nordgren, Phys. Rev. B 58, 1928 (1998). 
[34] T. H. Gfroerer, L. P. Priestley, F. E. Weindruch, and M. W. Wanless, Appl. Phys. Lett. 

80, 4570 (2002). 

[35] M. L. Benkhedir, M. S. Aida, A. Stesmans, and G. J. Adriaenssens, J. Optoelectron. Adv. 
Mater.7, 329 (2005). 

[36] A. Koizumi, H. Moriya, N. Watanabe, Y. Nonogaki, Y. Fujiwara, and Y. Takeda, Appl. 

Phys. Lett. 80, 1559 (2002). 
[37] H. Hirayama, A. Kinoshita, and A. Hirata, Y. Aoyagi, Appl. Phys. Lett. 80, 1589 (2002). 
[38] D. Penn, Phys. Rev. 128, 2093 (1962). 

[39] L. J. Wang, A. Kuzmich and A. Dogariu, Nature 406, 277 (2000). 
[40] D. Mognai, A. Ranfagni, and R. Ruggeri, Phys. Rev. Lett. Phys. Rev. Lett. 84, 4830 
(2000). 

[41] A. H. Reshak, Z. Charifi, and H. Baaziz, Eur. Phys. J. B 60, 463 (2007). 



14 

[42] A. Bhattacharyya, S. Lyer, E. Iliopoulos, A. V. Sampath, J. Cabalu, and I. Friel, J. Vac. 

Sci. Technol. B 20, 1229 (2002). 
[43] T. Someya, and Y. Arakawa, Appl. Phys. Lett. 73, 3653 (1998). 




O 0.25 0.5 0.75 

x(concentration) 



Fig.l Variation of direct bandgap with the concentration of x in Ali_ x Ga x N (0 < x < 1) 
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Fig.2 Frequency dependent dielectric functions of Ali. x Ga x N (0 < x < 1) 
(Real part: black lines, imaginary part: magenta lines) 
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Fig.3 Zero frequency limit of dielectric functions vs. direct bandgap of Ali. x Ga x N (0 < x < 1) 
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Fig.4 Frequency dependent refractive indices of Ali_ x Ga x N (Real part: black lines, imaginary 

part: magenta lines) 



o 

5 



1.9 



1.7 



3.5 



4.5 



Fig.5 Variation of static refractive index with bandgap of Ali_ x Ga x N (0 < x < 1) 
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Fig.6 Frequency dependent reflectivity of Ali_ x Ga x N (0 < x < 1) (Real part: black lines, 

Imaginary part: magenta lines) 



